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SUMMARY 


Through  the  use  of  linear  programming  techniques,  the  optimal  number  and  positioning  of 
patient  evacuation  assets  within  a  theater  of  operations  may  be  determined  to  ensure  the  orderly 
transport  of  casualties  from  the  front  lines  to  third  echelon  medical  treatment  facilities.  The 
Probabilistic  Location  Set  Covering  Problem  has  been  chosen  as  the  core  module  for  a  linear 
programming  model  to  assist  in  these  determinations.  The  Optimal  Placement  of  Casualty  Evacuation 
Assets  (OPTEVAC)  model  prompts  the  user  to  enter  the  dimensions  of  the  theater,  troop  deployment 
nodes,  types  of  evacuation  assets  available,  and  preferred  locations  of  medical  treatment  facilities. 
The  OPTEVAC  model  then  provides  output  as  to  the  required  numbers  of  ground  and  air  ambulances 
as  well  as  the  optimal  positioning  of  those  evacuation  assets  and  ambulance  exchange  points. 
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The  Optimal  Placement  of  Casualty  Evacuation  Assets: 
A  Linear  Programming  Model 


Combat  casualties  in  a  theater  of  operations  are  treated  at  mobile  medical  facilities  organized 
into  a  series  of  echelons  with  the  facilites  at  the  forward  echelons  having  the  greatest  mobility  but 
least  surgical  capability^  The  efficient  evacuation  of  casualties  from  echelon  to  echelon  is  essential 
to  ensure  the  wounded  personnel  reach  a  facility  with  the  capability  to  render  the  required  level  of 
care.  Distances  between  medical  treatment  facilities,  as  well  as  factors  such  as  the  type  of  terrain  and 
mode  of  transportation,  may  all  impact  the  evacuation  process.  Likewise,  the  number  of  treatment 
facilities  deployed  and  where  they  are  located  greatly  affects  the  casualty  handling  process  and  the 
ability  to  provide  adequate  casualty  care. 

In  the  Kuwaiti  theater  of  operations,  treatment  facilities  were  spaced  further  apart  than  the 
distances  considered  under  the  existing  medical  evacuation  doctrine;  consequently,  some  casualties 
requiring  evacuation  were  out  of  helicopter  range  or  in  locations  not  accessible  to  fixed-wing 
aircraft^.  While  casualties  were  relatively  few  in  the  Gulf  Conflict,  and  the  wounded  personnel  did 
receive  needed  treatment  in  a  timely  fashion,  future  operations  may  yield  larger  numbers  of  casualties 
with  a  greater  potential  for  overwhelming  the  casualty  evacuation  system.  Accurate  assessments  of 
the  required  numbers,  types,  and  deployment  locations  of  the  evacuation  assets  ensure  that  wounded 
personnel  are  transported  expeditiously  from  the  point  of  injury  through  the  various  echelons  of  care 
to  an  Echelon  III  treatment  facility,  as  dictated  by  the  severity  of  the  wound. 

The  focus  of  this  report  is  the  selection  and  development  of  an  appropriate  linear 
programming  technique  which  will  allow  optimization  of  the  number  and  positioning  of  evacuation 
assets  at  given  treatment  facilities  in  a  theater  of  operations.  Because  inter-theater  evacuation  of 
patients  is  handled  for  all  services  by  the  U.S.  Transportation  Command  (TRANSCOM),  the  present 
effort  focuses  on  intra-theater  casualty  evacuation  —  transporting  wounded  Marines  from  the  point 
of  injury,  through  Echelons  I  and  11  in  the  combat  zone,  and  ending  at  an  Echelon  III  level  of  care 
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(e.g.,  Fleet  Hospital).  The  methodology  developed  will  be  flexible  enough  to  allow  incorporation  of 
treatment  facilities  and  evacuation  assets  of  other  service  branches  within  the  combat  zone, 
particularly  those  of  the  Army  which  may  very  well  have  forces  interdeployed  with  the  Marines. 

The  optimization  planning  tool  will  incorporate  the  projection  methodology  of  the  FORECAS 
system^  which  simulates  daily  casualty  averages  and  maximum  daily  patient  loads.  The  simulations 
generated  by  FORECAS  are  based  on  the  statistical  properties  of  empirical  data  observed  from  four 
previous  ground  combat  operations'*.  In  addition  to  patient  flows,  the  optimization  model  will  require 
data  on  troop  strengths  and  deployment  locations,  terrain,  size  of  theater,  numbers  and  types  of 
ground/air  ambulances  available,  types  and  locations  of  Echelon  II  and  Echelon  III  treatment 
facilities,  and  length  of  evacuation  delay.  Using  this  input,  the  planning  tool  will  compute  the 
optimum  numbers,  types,  and  deployment  coordinates  of  the  evacuation  assets,  including  ambulance 
exchange  points. 

A  comprehensive  examination  of  existing  optimization  models  was  conducted  and  the 
Probabilistic  Location  Set  Covering  Problem  (PLSCP)’’^  was  determined  to  be  the  one  that  most 
accurately  represented  the  demands  and  constraints  needed  to  model  the  intra-theater  evacuation  of 
casualties.  Other  optimization  models  that  were  examined  include  the  location  set  covering  problem 
(LSCP),  the  maximal  covering  location  problem  (MCLP),  the  maximum  availability  location  problem 
(MALP),  and  two  derivatives  of  the  chosen  PLSCP:  the  a-reliable  p-center  problem  and  the 
maximum  reliability  location  problem  (MRLP).  The  PLSCP  was  chosen  over  the  other  models 
because  it  minimizes  the  total  number  of  required  evacuation  assets  within  the  theater,  while  at  the 
same  time  incorporating  randomness  in  vehicle  availability.  The  vehicles  of  the  other  models 
investigated,  are  sited  only  on  the  basis  of  geographical  coverage  and  not  on  the  basis  of  availability. 
The  assumption  underlying  these  alternative  models  is  that  the  vehicles  are  continuously  available  for 
evacuation  of  casualties.  In  low  demand  operations,  this  assumption  is  not  unreasonable;  however 
in  high  tempo  military  operations  where  frequent  casualty  calls  keep  evacuation  assets  on  the  roads, 
this  assumption  cannot  be  Justified.  Use  of  the  PLSCP  algorithms  within  the  evacuation  model 
(OPTEVAC)  will  allow  medical  planners  to  input  the  “alpha-level”  at  which  the  optimization  model 
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is  to  mn  ~  a  statistic  which  represents  the  percentage  of  the  maximum  daily  casualty  load  for  which 
the  proposed  assets  will  meet  the  demand.  The  maximum  daily  casualty  load  will  be  computed  by 
adding  the  maximum  projected  wounded-in-action  (WIA)  incidence  to  the  maximum  projected 
disease  and  non-battle  injury  (DNBI)  incidence  during  the  operation.  The  alpha  level  is  important 
because  of  the  casualty  “spikes”  that  occur  during  a  combat  operation.  For  instance,  the  planner  may 
want  to  economize  resources  by  deploying  dedicated  assets  to  ensure  evacuation  coverage  for  80% 
of  the  maximum  casualty  load  (a  =  .80)  and  rely  on  “vehicles  of  opportunity”  to  transport  the 
casualties  that  exceed  the  capabilities  of  these  dedicated  assets.  Alternatively,  if  evacuation  assets 
were  programmed  for  the  maximal  projected  daily  casualty  load,  many  assets  would  be  used  but  a 
small  percentage  of  the  time,  if  at  all. 

Patient  flows  within  the  OPTEVAC  model  will  be  determined  from  the  user-specified  input 
of  troop  placement,  with  projections  based  on  casualty  distribution  trends  evidenced  during  previous 
combat  operations.  The  patient  stream  is  actually  composed  of  two  separate  medical  admission 
flows:  WIA  and  DNBI.  The  daily  WIA  flow  is  generated  by  drawing  a  random  deviate  from  an 
exponential  distribution  based  on  ‘battle  intensity’ -specific  averages  and  geographical  theater 
considerations.  Projected  DNBI  incidence  is  similarly  derived  using  variates  drawn  from  a  lognormal 
distribution.  The  casualty  generator  also  incorporates  two  other  statistical  patterns  observed  within 
the  empirical  data:  cross-correlation  betwen  WIA  and  DNBI  and  autocorrelation  within  the  WIA 
flows^^. 

The  following  section  describes  the  linear  programming  model  (PLSCP)  chosen  for  optimizing 
casualty  evacuation  assets  and  provides  a  detailed  analysis  of  how  the  problem  is  mathematically 
represented.  Essential  to  accurate  assessment  of  the  numbers  of  evacuation  assets  required  are 
variables  relating  to  the  specific  types  of  assets  (transportation  mode,  litter  capacity,  range,  mobility, 
etc.)  as  well  as  factors  related  to  terrain  and  weather.  The  scope  of  the  present  report,  however,  is 
the  linear  program  to  be  used  in  optimizing  evacuation  assets;  the  parameters  relating  to 
environmental  conditions  and  ambulance  characteristics  are  the  subject  of  a  parallel  effort. 


5 


THE  MODEL 


The  OPTEVAC  model  will  provide  medical  planners  with  the  required  number  and  optimal 
placements  of  evacuation  assets  to  ensure  sufficient  casualty  transport  while  minimizing  oversupply 
of  ground  and  air  ambulances.  The  Probabilistic  Location  Set  Covering  Problem  (PLSCP)  linear 
program  has  been  selected  as  the  core  module  for  OPTEVAC,  and  it  has  been  modified  to  make  it 
applicable  to  the  optimization  of  casualty  evacuation  assets  in  a  combat  theater. 


The  PLSCP  module,  incorporated  into  the  casualty  evacuation  application,  seeks  the  positions 
of  z  evacuation  assets  in  a  theater  of  operations  which  minimizes  the  number  of  vehicles  (ground  and 
air  ambulances)  required  so  that  there  is  a  95%  certainty  that  the  user-specified  evacuation  demands 
will  be  met.  Mathematical  formulation  of  the  PLSCP  for  the  casualty  handling  process  is; 
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H  = 

Variables; 
z  = 

a  = 

j,j  = 

dji  = 

S 

Nj 

bi 


frequency  of  casualty  calls,  or  trips,  at  demand  node  k  (casualty  calls/day); 
the  set  of  demand  nodes  within  S  of  demand  node  i. 

the  total  number  of  evacuation  assets  distributed  over  all  of  the  facilities  in 

the  theater  of  operations  by, 

index  and  set  of  demand  nodes; 

index  and  set  of  treatment  facility  sites; 

the  distance  from  facility  site  j  to  demand  node  i; 

the  distance  standard  within  which  a  treatment  facility  is  desired  to  be  found; 
{jjdji  <  S};  the  set  of  facility  sites  within  S  of  demand  node  i; 
the  number  of  evacuation  assets  at  facility  site  j; 
the  minimum  number  of  evacuation  assets  required  within  M; . 


The  OPTEVAC  planning  tool  consists  of  input  screens  which  employ  graphical  user  interfaces 
prompting  the  user  for  the  information  needed  for  the  simulations.  The  first  input  that  the  planner 
must  supply  into  the  OPTEVAC  model  is  the  dimensions  of  the  theater.  An  image  of  a  grid  with  the 
user-specified  dimensions  will  then  be  displayed  on  the  computer  screen  immediately  following  this 
input.  The  next  required  input  pertains  to  location  of  troop  deplo5mient.  Following  troop 
deployment  input,  the  user  will  be  asked  to  specify  the  numbers  and  types  of  medical  treatment 
facilities  and  to  indicate  their  locations  within  the  theater.  Other  key  inputs  to  the  model  will  include 
factors  such  as  battle  intensity,  geographical  region,  expected  length  of  the  operation,  evacuation 
delay,  “alpha-level”  associated  with  the  maximum  projected  daily  casualty  load,  and  types  of 
evacuation  assets  available. 


The  standard  PLSCP  linear  program,  by  itself,  is  not  sufficiently  robust  to  handle  the  demands 
of  the  OPTEVAC  model,  so  modifications  were  required  to  fit  it  to  the  casualty  evacuation  model. 
The  PLSCP  was  originally  written  to  satisfy  the  demands  of  a  city  population  in  need  of  ground 
ambulances  in  a  timely  fashion.  The  urban  transport  problem  assumed  only  one  type  of  ambulance, 
which  implied  one  speed  and  one  litter  capacity.  In  a  military  theater,  however,  there  is  a  multitude 
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of  vehicles  available,  and  each  has  it’s  own  distinct  speed,  range,  and  litter  capacity.  Furthermore, 
within  a  combat  theater,  different  types  of  vehicles  need  to  be  dispatched  simultaneously. 

To  integrate  more  than  one  type  of  evacuation  asset  into  the  model,  information  regarding 
the  capacities,  speeds,  and  the  numbers  of  each  vehicle  available  for  deployment  must  be  included. 
Both  ground  and  air  ambulances  are  actively  involved  in  the  casualty  evacuation  process,  so  a  revised 
algorithm  is  needed  to  determine  the  demand,  or  frequency  of  calls  per  day,  of  a  given  configuration 
of  transportation  assets.  This  algorithm,  which  represents  the  core  of  the  OPTEVAC  model, 
considers  the  placement  and  allocation  of  evacuation  assets  for  varying  litter  capacities. 

The  model  performs  the  necessary  distance  calculations  between  the  troop  deployment  nodes 
and  the  medical  treatment  facilites  and  begins  to  place  the  various  types  of  available  vehicles  at 
treatment  facilites  to  solve  for  the  optimal  vehicle  configuration.  Different  placements  of  troop  nodes 
and  treatment  facilites  will  of  course  yield  different  numbers  and  combinations  of  each  vehicle  located 
at  the  facilities.  Not  always  will  all  of  the  available  vehicles  be  required  for  a  given  treatment  facility 
configuration  since  the  model  finds  the  minimum  number  of  required  vehicles  to  handle  the  casualty 
flow.  Similarly,  the  results  may  indicate  that  more  vehicles  must  be  made  available  than  were 
originally  entered  to  handle  the  casualty  load. 

The  OPTEVAC  program  utilizes  the  PLSCP  module  to  solve  for  the  optimal  number  of 
evacuation  assets  given  a  specific  configuration  of  treatment  facilities.  As  part  of  the  PLSCP  model, 
the  first  task  is  to  determine  each  of  the  b;  variables  (minimum  number  of  required  assets  within 
demand  area),  and  there  will  be  as  many  b;  variables  as  there  are  demand  nodes.  Each  b,  can  be 
solved  after  the  planner  has  entered  the  aforementioned  required  inputs  into  the  program,  and  solving 
for  each  b;  will  be  accomplished  by  using  Newton’s  Method*.  Newton’s  Method  will  be  used  to  find 
the  bj’s  after  the  user  has  entered  the  inputs  required  to  solve  for  F;  (total  evacuation  asset  demand). 
The  required  bj’s  will  be  inserted  into  the  constraint  matrix  of  the  PLSCP  and  hence  complete  the 
system  of  equations.  There  will  be  “i”  (#  of  demand  nodes)  equations  and  “j”  (#  of  treatment 
facilities)  variables  in  the  constraint  matrix. 


8 


The  following  casualty  evacuation  problem,  solved  by  the  Simplex  Method®,  demonstrates 
that  OPTEVAC  can  yield  solutions  for  multiple  casualty  evacuation  asset  types.  In  this  facile 
example,  consider  a  60-day  operation  in  which  there  are  five  Echelon  II  facilities,  three  troop 
deployment  nodes  with  strengths  seen  in  Figure  1,  and  that  the  average  WIA  and  DNBI  rates  are  5.77 
per  1000  strength  per  day  and  4.22  per  1000  strength  per  day,  respectively.  Also  assume  that  10 
ground  ambulances  and  2  air  ambulances  are  available  for  deployment.  The  first  calculation 
performed  is  to  determine  the  frequency  of  calls  per  day  at  each  troop  deplojnnent  node  within  the 
theater.  The  necessary  algorithm  must  factor  in  the  placement  and  allocation  of  evacuation  assets  for 
varying  litter  capacities.  Therefore,  the  PLSCP  algorithm  needs  to  be  modified  such  that  separate 
fi-equencies  (^ground  and  and  separate  demand  areas’  vehicles  (bi^grou„d  and  are  computed 
for  each  type  of  vehicle.  Defined  below  is  an  application  of  the  algorithm  that  calculates  the 
frequency  of  calls  per  day,  when  multiple  vehicle  assets  are  being  used  for  evacuation  (vehicle 
specifications  are  also  listed).  For  instance,  with  only  two  types  of  vehicles  (e.g.,  a  single  ground 
ambulance  type  and  a  single  air  ambulance  type)  this  algorithm  computes  fk^ground  and  fk,airas  follows: 


Litter  capacity 

Speed* 

Quantity  available 

Ground  ambulance 

4 

16  km/hr 

10 

Air  ambulance 

6 

150  km/hr 

2 

*  The  speeds  of  the  dedicated  assets  include  the  loading  and  unloading  of  patients*®. 
Define  variables: 


a/g  ratio 


air  ambulance  speed 
ground  ambulance  speed 


150  km! hr 
16  km! hr 


9.375  , 


g  capac 
a  capac 
g  avail 
a  avail 
tot  avail 


=  ground  ambulance  capacity, 

=  air  ambulance  capacity, 

=  number  of  ground  ambulances  available, 

=  number  of  air  ambulances  available, 

=  the  total  number  of  ground  and  air  ambulances  available. 
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g  factor 

2  avail  , 

^  (S  capac)  , 

tot  avail 

a  factor 

a  avail  ,  x  ^  , 

-  X  {a  capac)  x  (a/g  ratio) 

tot  avail 

Z  factor 

_  maximum  daily  casualties 

g  avail  +  a  avail 

'‘■k.ground 

^  (Z  factor)  X  {g  factor) 

g  capac 

4air 

_  (Z  factor)  X  {a  factor) 

? 

a  capac 


The  sum  of  the  frequency  of  calls  from  all  demand  nodes  within  the  theater  (fk,grou„<,  +  4^)  will  be  used 
to  calculate  F;,  which  will  then  lead  to  the  calculations  of  the  b/s.  By  way  of  the  Simplex  Method, 
the  PLSCP  module  determines  that  3  ground  ambulances  and  1  air  ambulance  are  required  within  Mj 
(demand  area  1),  4  ground  ambulances  and  1  air  ambulance  are  required  within  M,  (demand  area  2), 
and  4  ground  ambulances  and  1  air  ambulance  are  required  within  Mj  (demand  area  3)  as  well 
(symbolically,  for  ground  ambulances:  =  3,  =  4,  4, ground  =  4,  and  for  air  ambulances: 

^l.air  ^2, air  ^3, air  0- 

Referring  to  Figure  1,  the  corresponding  mathematical  setup  of  the  PLSCP,  when  solving  for 
the  optimal  number  of  ground  ambulances  in  the  example  scenario,  is  the  following: 
minimize  z  =  Xu  +  X21  +  Xj,  +  X41  +  X51 
subject  to 

^11  ^21  -  3  (troop  node  1  is  within  range  of  MTF  1  and  2) 

^21  ^31  ^41  ^51  ^  4  (troop  node  2  is  within  range  of  MTF  2,  3,  4,  and  5) 

X4,  +  Xjj  >  4  (troop  node  3  is  within  range  of  MTF  4  and  5) 

where  Xj  are  positive  integers, 

and  when  solving  for  the  optimal  number  of  air  ambulances  in  the  example  scenario,  the 
corresponding  mathematical  setup  is  the  following: 


10 


Kilometers 


Figure  1 .  Sample  input  screen  to  the  Optimization  of  Casualty  Evacuation  Assets 
(OPTEVAC)  model  showing  troop  deployment  nodes  (X's)  and  medical 
treatment  facilities  (+'s). 

1 1 


minimize  z  =  Xj,  +  X22  +  Xj,  +  X42  +  X52 

subject  to 

Xj2  +  X22  +  X32  +  X42  +  Xj2  ^  1  (troop  node  1  is  within  range  of  each  MTF) 

X12  +  X22  +  X32  +  X42  +  ^^52  -  ^  (troop  node  2  is  within  range  of  each  MTF) 

X|2  +  X22  +  X32  +  X42  +  Xj2  ^  1  (troop  node  3  is  wdthin  range  of  each  MTF) 

where  Xj  are  positive  integers. 

In  the  aforementioned  equations,  the  first  subscript  represents  the  Echelon  II  treatment 
facilities,  and  the  second  represents  the  type  of  vehicle.  The  first  equation  of  the  ground  ambulance 
constraint  matrix  may  be  interpreted  as  “troop  deployment  node  1  is  within  range  of  medical 
treatment  facilities  1  and  2  by  way  of  ground  ambulance”.  Similarly,  each  equation  in  the  air 
ambulance  constraint  matrix  indicates  that  none  of  the  troop  deployment  nodes  are  out  of  air 
evacuation  range  for  each  of  the  five  Echelon  II  treatment  facilities. 

After  solving  the  above  linear  programming  problems,  the  optimal  solution  was  found  to  be 
placement  of  two  ground  ambulances  at  treatment  facility  1,  one  at  treatment  facility  2,  one  at 
treatment  facility  3,  three  at  treatment  facility  4,  and  one  at  treatment  facility  5,  for  a  total  of  eight 
necessaiy  ground  ambulances  (x„  =  2,  X2,  =1,  X31  =  1,  X4,  =  3,  X5,  =  1,  and  =  8).  By  inserting 
these  results  into  the  constraint  matrix  of  the  PLSCP  for  ground  ambulances,  it  is  clear  that  the  results 
are  indeed  a  minimum; 

2  +  1  >3 

1  +  1  +  3  +  1  >  4  => 

3  +  1  >  4 

thus  the  constraint  equations  hold,  and 

Aground  ~  ^11 ''’^21 '*’^31 '*'^41  "'’^51  ~  2  +  1+1 +3+1  =  8. 

Further,  only  a  single  air  ambulance  was  required  (z^  =  1)  and  could  be  placed  at  any  Echelon  II 
facility,  though  placing  it  in  the  center  facility  would  be  the  logical  choice. 


3  >  3 
6  >  4 

4  >  4 
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In  this  sample  scenario,  the  optimization  algorithms  were  able  to  minimize  the  number  of 
evacuation  assets  required  in  the  operation.  Ten  ground  ambulances  and  two  air  ambulances  were 
available  for  the  operation,  but  only  eight  land  transports  and  one  air  asset  were  necessary  to  perform 
the  evacuation  demands  of  the  defined  scenario.  OPTEVAC  finds  the  optimum  numbers  of  required 
dedicated  assets  for  a  given  theater  of  operations,  and  in  the  example  scenario,  the  planner  saved  on 
resources  by  20%  for  ground  ambulances  and  50%  for  air  ambulances.  These  results  demonstrate 
that  the  OPTEVAC  model  allows  the  user  to  determine  the  minimum  numbers  and  optimal 
positioning  of  patient  transportation  assets  required  to  meet  the  casualty  handling  demands  of  the 
operation. 


Summary 

The  PLSCP  model  is  the  appropriate  choice  for  modeling  the  efficient  transport  of  casualties 
between  intra-theater  medical  facilities.  The  model  will  be  integrated  into  a  software  environment 
that  will  allow  planners  to  project  the  optimal  numbers  and  placement  of  assets  needed  for  casualty 
handling  in  a  theater  of  operations.  The  user  will  be  able  to  input  the  dimensions  of  the  theater,  the 
quantity,  types,  and  placement  of  medical  treatment  facilities,  as  well  as  the  types  of  evacuation  assets 
available.  The  OPTEVAC  casualty  evacuation  model  will  then  provide  the  optimal  number  and 
position  of  the  needed  ground  and  air  ambulances. 
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